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There are two related theorems which hold even in far from equilibrium, namely fluctuation 
theorem and Jarzynski equality. Fluctuation theorem states the existence of symmetry of fluctuation 
of entropy production, while Jarzynski equality enables us to estimate the free energy change between 
two states by using irreversible processes. On the other hand, relationship between these theorems 
was investigated by CrooksQ] for the classical stochastic systems. In this letter, we derive quantum 
analogues of fluctuation theorem and Jarzynski equality in terms of microscopic reversibility. In 
other words, the quantum analogue of the work by Crooks[l| is presented. Also, for the quasiclassical 
Langevin system, microscopically reversible condition is confirmed. 

PACS numbers: 05.70.Ln,05.40.-a,05.30.-d 



I. INTRODUCTION 

Since the discovery by Evans et al.0, many types of fluctuation theorem has been presented both for 
deterministic E l El [1 and stochastic systems Q [3 0] G3 ■ Experimental verification of fluctuation theorem was 
also performed[ll|. Though the systems, the definition of entropy production, and interpretations (such as distinction 
of transient or steady state) are different, fluctuation theorems have a following universal form. : 

P(AS = A) m 



P(AS = -A) 



Here, AS is the entropy generated, and P(AS = A) is the probability for AS — A. A few quantum analogues of 
fluctuation theorem has been proposed ^3 H3 [13 • I n the probability for the energy change between initial and 
final states is considered, so we need to perform observations twice, while in |l l| the spectrum distribution of work 
operator is considered and we need to perform only one observation. In |l5j| . probability for heat exchanged between 
two systems is considered, and we need to perform observations twice. The features of the two times measurement 
approach and the work operator approach are discussed in |2l| . 

In 1997, Jarzynski presented nonequilibrium equality for free energy difference(Jarzynski equality) [l6j | which is 
a generalization of the minimum work principle. Jarzynski equality states that the change of free energy AF is 
calculated by taking the sample average of exponentiated work done on a system e~@ w along irreversible processes 
which has a following form: 

AF = -ibg(e-^> (2) 

Here, AF is the free energy difference, (3 is the inverse temperature, W is the work done on a system along some 
irreversible process, and () denotes the sample average. The remarkable feature of this equality is that the change 
of equilibrium quantity can be estimated by using the irreversible process. Jarzynski equality was experiment ally 
confirmed by RNA stretching experiment 17]. Quantum analogues of Jarzynski equality was proposed by Mukamel|l8j] 
and by Yukawa 19]. Both quantum analogues employ the situation where a system is initially in thermal equilibrium 
and obeys the canonical ensemble, and the time evolution is described by the time dependent Hamiltonian H(t). In 
[l8l | Jarzyski equality is derived in terms of master equation, while in 0] the density matrix approach is considered. 
Here we propose a derivation based on the microscopic reversibility (JSJ. 

For classical stochastic systems, Crooks Jj derived classical fluctuation theorem and Jarzynski equality in terms of 
microscopic reversibility. In this letter, we give a unified derivation of the quantum analogue of fluctuation theorem 
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and Jarzynski equality along the thought of Crooks [lj. The ratio of the probability distributions of entropy generated 
between time forward and time reversal process is considered, and we define the entropy production analogous to that 
of 0. So, fluctuation theorem derived here can be considered as a quantum extension of that of pj. 

Our emphasis is on the unified treatment of these theorems in terms of microscopic reversibility. This condition is 
derived for a quasiclassical Langevin system in the appendix. 



II. QUANTUM FLUCTUATION THEOREM 
A. the case of thermally isolated system 

In this subsection, we derive a quantum analogue of fluctuation theorem for a system which couples to an external 
driving, thermally isolated and whose Hamiltonian H (t) is invariant under time reversal. In order to perform this, we 
derive the relation (J3J). In 0] for the case of time independent system, the relation © was discussed, however, this 
relation is necessary for the derivation of fluctuation theorem here, we show the derivation. 

Let the unitary time evolution operator between t = and t = T be Ut- And let the density matrix at t = — 
be p(0) and define p(T) = U T p(0)U T - Consider some arbitrary observables A and B whose eigenstates \a n ) and \b m ) 
form the normalized orthogonal complete set. 

Then, we perform the measurements about A at t = +0 and about B at t — T. Firstly, we show the following 
relation. 

-Pf(K) -> \b m )) = log<a„|p(0)| a „>-log<6 m | / >(T)|6 m ) foN 

p R (e\b m ) - eK» ' [ 1 

Here, Pf {\dn) ^ \b m )) is the probability that the forward process II whose observed states at time t = and t = T 
are \a n ) and \b m ) occurs. O is the antilinear time reversal operator, and Pii(@\b m ) — ► 0|a n )) is the probability that 
the time reversed process II* occurs, namely at t = Q, we perform the measurement about B and find the initial 
state as 0|6 m ) according to the ensemble with the time reversed density matrix 0p(T)O and at t = T perform the 
measurement about A and find the final state as 0|a„). We note that in ^| the average over PfQolq) — > \ar}) of 
the logarithm of the ratio p^Q\aT)-^&\al)) ^ s sri0wn to be equal to the microcanonical entropy change in the case |a ) 
and |ar) are macro states provided the Hamiltonian is time independent. 

The derivation of 10 is as follows. Because the probability we find the initial state as \a n ) is (a n \p(0)\a n ), one has 

Pf(K) -> \b m )) 



P R (0\b m ) -> S\a n )) 

(a»|p(0)|a»)|(6 m |[/ T |a»)| 2 

(fe m |eep(T)ee|6 m )|(a„|e^e|6 m })| 2 

_ e log<a„|p(0)|a„>-log(b m |p(T)|b m >) 

Here, as mentioned above, Ut is the unitary time evolution operator and its time reversal is U T , express that anti 

linear operator acts to the left, and we used the relation for unitary operator (b m \UT\a n ) = {a n \ ®U T Q\b m ). This 
relation is proved as follows : Because we assumed that the Hamiltonian is time reversal invariant, QH(t) — H(t)Q, 
one has 

(a n \QU T Q\b m ) 

= lim (a n \QU\U2--U Ar0|m) 

= (b m \U T \a n ) (5) 

Here, Uk denotes the unitary time evolution operator between time ^j^-T and j?T. 

Then one has the quantum analogue of fluctuation theorem for the system without heat bath as shown below. 
We define a counterpart of the entropy production AS analogous to the classical entropy production defined in pj. 

AS = fc B log(a„|p(0)|a„) - fc B log(6 m |p(T)|6 m ) (6) 

Of course this definition of entropy production depends on the choice of the observables A and B. The state is projected 
into other basis by measurement and different ways of measurements may cause different entropy productions, unlike 
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in classical system where all observables commute with each other and the value of the entropy production is unique. 
We note that AS is clearly considered as entropy production when these observables diagonalize the density matrices 
p(0) and p(T) respectively. For example, when the system is in equilibrium at initial and final time and we choose A 
and B as Hamiltonian, this condition is satisfied. We treat this case in the section V. In such cases, AS is considered 
as entropy production in the sense of namely AS = — k B log pf + ks logpi — Here, pi/f is probability density 
of initial/final times, Q is the heat transferred into the system from the heat bath, Tg is the temperature of the heat 
bath. In comparison to the diagonal representation of von Neumann entropy S = — Tr(plogp) = — £ n pn logp n , the 
quantity — ks log/9/ + ks logpi is considered as entropy production of the system, and — is the entropy production 
of the heat bath. From the assumption that the system is thermally isolated, heat transferred to the system Q is 0. 
Then, one obtains a quantum analogue of fluctuation theorem : 

Pf(AS) = E n>m P F (\a n ) -> \b m ))d(AS - fc B (log(a„|p(0)|a n ) - log(b m \p(T)\b m })) 

= e^i^ m P R {Q\b m ) -> Q\a n ))S(AS + k B {\og(b m \p(T)\b m ) - log<a n |/>(0)|a„))) 

= P R (-AS)e$i. (7) 

Here, P F (AS)(P R (AS)) is the probability that the counterpart of entropy production along the forward(backward) 
process is AS". 



B. the system with heat bath 

In this subsection, we consider a system coupled to a heat bath whose total Hamiltonian is H (t) = H s (t) + H R + Hi- 
Here, H s (t) 7 Hb, and Hi are the Hamiltonian of system, heat bath and interaction respectively. Through out this 
subsection, we choose A and B as system Hamiltonian H s (0) and H S {T) respectively. And we denote the n'th 
eigenvalue of H s (t) as E n (t) and corresponding eigenvector as \n,t). We restrict ourselves to the system where the 
microscopic reversible condition is satisfied. 

P F (\n,0) ^ \m,T)) _ fiQ 



P R (Q\m,T) ->e|n,0)) 



e-^ (8) 



Here, P F (\n, 0} — > |m, T)) denotes the conditional probability that the system is initially in the energy level \n, 0) and 
after some duration time T jumps to the level \m,T). The distinction of time forward and reversed processes is the 
same as previous section. Q is the heat absorbed by the heat bath. This relation plays an essential role in for 
the derivation of Fluctuation Theorem. Also in quantum system, there are many situations where the microsocopic 
reversibility does hold. In appendix, as an physically important system which satisfies this relation, we confirm this 
microscopic reversibility for the so-called quasiclassical Langevin system |2j| . The rest of this section is devoted for the 
derivation of fluctuation theorem for the system with the heat bath in terms of the relation JSJ. The total probability 
Pp ot (\n, 0} — ► 1 77i, T)) that the system is initially in \n,0) and after some duration T found in |m, T) is given as the 
product of initial state probability and conditional probability. 



?Toti^ , ,, ,, , ,,,,,„ ,, a* 

e PQ = e k B 



P£ ot (\n, 0} \m, T» _ (n, 0\p a (0)\n, 0) 



P| ot (9|m,r) — 6|n,0)) (m,T\eep s (T)eQ\m,T) 
AS 

— = log(n, 0|p s (0)|n, 0} - log(m, T\p s (T)\m, T) ~ f3Q. (9) 
k B 

Here, Q is the heat transferred to the system from the heat bath, and defined the reduced density matrix p s (t) = 
TrB(pTot) as a system density matrix . 

Then one has quantum analogue of fluctuation theorem for the system with a heat bath 

AS 

P F (AS) = E„. m Pj ot (|n,0) -> \ m ,T))8(— - (log(n, 0|p s (0)|n, 0) - log(m, T\p s (T)\m, T) - 0Q)) 

kb 

A S rr 1 AS 

= e^Z n , m P£ ot (<d\m, T) e n, 0) )S (-=— + (log (m, T\p s (T) m, T) - log(n, \p s n, 0} + PQ)) 

Kb 

= P R (-AS)e^i. (10) 
This relation corresponds to the classical fluctuation theorem derived in 0. 
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III. QUANTUM JARZYNSKI EQUALITY 

In this section, we derive quantum Jarzynski equality in terms of microscopic reversibility JSJ . 

At first, we set up the framework on which we discuss here. We assume that the heat bath is large enough and 
interaction with the system is weak enough so that the following conditions do hold. As the system density matrix, 
we use the reduced density matrix p s (t) = TrB(pTot(t))- In order to discuss the free energy change between two 
equilibrium states, we require that the density matrices of the system p s (0) and p s (T) are canonical distribution at 
the same temperature. We set observables A and B as the system Hamiltonian H s (0) and H S (T) as previous section 
and denote the eigenstate of H s (0) and H S {T) as |n, 0} and |m,T) with the eigenvalues -EVi(O) and E m {T). Then 

(n,0|p a (0)|n,0> = -L e -/3(«>oi^(o)Ko> 
Z\ 

(m,T\p s (T)\m,T) = ± e -0MH.(T)\ m ,T) 

(11) 

Here, Z\ and Z^ are the partition functions for p s (0) and p s (T). Therefore, the relation © becomes as 

P^(\n,0)^\m,T)) 



pT ot (@\m,T) 6|n,0» 

^2 p(( m ,T\H,(T)\m,T}-(n,O\He(0)\n,0}-Q) 
_ e /3(E m (T)-E„(0)-O)-AF) / 12 s 

Here, AF = — i log ^ is the change of Helmholtz free energy. We note that W = E m {T) — E n (0) — Q is considered 
as work externally done on the system. So one has 



(13) 



Pj 0t (M)^|m,T)) = MW-AF) 

pT ot {0\m,T) -> e|n,0» 
Taking average of both sides, one has 

(e-^- AF) ) = WJ° 4 (M) - \m,T)) ^§^^^ 

Pp ot {\n,0) -> \m,T)) 

= 1. (14) 
Here, (..) denotes the average over the probability Pj°*(|n,0) -> \m,T)). Therefore, one has 

Air = _Ibg(e-^> (15) 

This is the relation that we call the quantum Jarzynski equality. 

In summary, we derived a quantum extension of Jarzynski equality and fluctuation theorem in terms of microscopic 
reversibility ©■ And this relation |JSJ is confirmed for quasiclassical Langevin system. This unified treatment is the 
quantum version of Crooks' derivation of fluctuation theorem and Jarzynski equality for classical system. 
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APPENDIX A: MICROSCOPIC REVERSIBILITY FOR QUASICLASSICAL LANGEVIN SYSTEM 

As an physically important system which satisfy the microscopic reversibility (JSJ, we confirm this relation for the 
so-called quasiclassical Langevin system [23] ■ This section is devoted for the derivation of relation JSJ for the system 
which obeys the quasiclassical Langevin equation. For the classical Langevin system, microscopic reversibility was 
derived by Narayan and Dhar • Here, we utilize their result for quantum Langevin system in semi classical regime. 
To do this, it is most convenient to use the quantum noise theory by Gardiner |24j| . Because a straightforward 
derivation of Langevin equation is known |24j . we consider a system which interacts with a harmonic reservoir. Let 
the total Hamiltonian be 

H{t) = + V(q, A(t)) + ~ J d\«p x - n x qf + <Jtf x ) (Al) 

Here, q,p are the canonical coordinates of system and q\,p\ are those of heat bath. X(t) in the potential term is the 
control parameter corresponding to the external agent. Note that for simplicity, we consider the coupling between 
position and momentum, this is equivalent to the position-position coupling, which might be rather natural form. 
The system is assumed to be initially uncorrelated to the heat bath and the initial density operator of the heat bath 
Pb is assumed to be canonical. (p(0) = p s ® pb, Pb = -^e~P f dXhuJxt - a * ax + ^). This model Hamiltonian is standard 
and one of the ideal examples which describes the system interacting with the heat bath. Let Y(i) be an arbitrary 
system operator in Heisenberg picture. Firstly, we define the quantity p(t) by 

Tr s (Y(t)p s ® p B ) = Tr s (Y <g> p{t))p B (A2) 
The equation of motion for p(t) (adjoint equation) is given as 

(i(t) = A a p{t) + a(t)A 1 p{t) (A3) 
where A p(t) = j-[H sl p{t)} + i[[ 7 g, p{t)} + : q] , A 1 p(t) = £[?,ji(t)] and a{t)p{t) = i[£(t), M (t)]+. = 
i J d\n\^J ^Y^(—ax{0)e^' luJxt + a i x (Q)e luJxt ) is the Langevin force and here we assume the density of the state is con- 
stant — so that the adjoint equation be Markov (The Markov approximation). Then £(i) satisfies the Fluctuation 
Dissipation Theorem, 

([£(*)> f (01+) = — / dwwcothCSfiw) cos «(<-«'), (A4) 
^ Jo 

where the average (...) = Tre(...pB) denotes the average over the bath variables. This adjoint equation is rewritten in 
the form of Kramers equation. Suppose the Wigner function corresponding to p,(t) is W(q,p, t) = J dr(q + ^\p(t)\q — 
!)e -l T? . Then one has 

- dT = t- Tx i + ^ A(t)) + 7 £ - a(t)))W + ^-r—V*+H*, Ht)W (A5) 

Due to the associative nature of a(t), the first term of this equation is equivalent to the following c-number quasi- 
classical Langevin equation. Note that up to 0(h) the second term is negligible and further more, this approximation 
becomes better in the large friction limit. 

P 

x = — 
m 

p = -V'(x,\) - 7 H + fi(t) 
m 

Here a(t) is the c-number stochastic process whose all the moments are equal to that of a(t) respect to the thermal 
average. 

(a(t)a(t')) = Tr B (a(t)a(t')pB) etc. (A6) 

In particular, under the definition of thermal noise operator £(t) and the assumption that pb is canonical, we can 
show that a(t) is Gaussian process. Thus the mean and variance are sufficient to determine the stochastic feature of 
5(t). 

(«(*)) = o 

(a(t)a(t')) = — / dtuuj coth(PHuj) cos uj(t - t') 
n Jo 
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We derived qausiclassical Langevin equation, and then we can use the result of Narayan and Dhar 7] to obtain the 
relation (JSJ. From the Gaussian nature of a(t), one has the probability functional Pp/ft[a(t)] for each realization of 
a(t) 

Ppi R [5l(t)] = Ce~^fo dS f° ds ' & ^^ /o°° dwwcath( / 8/to)caB«(s-«'))- I S(s / ) 
= Ce~tti ds J ds'a(s)(^S(s~s')+0(h 2 ))a(s') 

A(s,s')~ 1 denote the inverse of A(s,s') as a kernel. Here the distinction of time forward and reverse F/R means 
that the time reversal of the control parameter X(t) is also considered. As mentioned above, we consider up to 0(H) 
(semi-classical regime) and thus the 0(h 2 ) term is omitted. Then after the same discussion of 0, one has the following 
relation. 

PF(qf,Pf\qi,Pi) = e -(3(H B (q f ,p f Mt))-H s ( q i,P.M0))-W) / A7 ) 

Pr{<H, -PiWf, ~Pf) 

Here, PF/R{qfTPf\qi,Pi) is the probability that the initial and final states for the quasiclassical Langevin equation 

2 

are (qi,Pi) and (qf,Pf) respectively. H s (q,p,\(t)) = f— + V(q,X(t)) is the system Hamiltonian at time t and 

W — J Q ds-^V(q, X(s))X(s) is the work externally done through the control parameter X(t). Then we choose the 
initial and final points of the phase space according to the initial and final Wigner functions. And the probability 
that the initial and final states are \n, 0} and m, T) is calculated as 

P F (\n,0) -> \m,T)) = ( J dqidpiW n (qi,Pi,0) J dq f dp f W m (q f ,p f ,T)S{q f - q{t 7 q l ,p l ))5(p f - p(t, q t ,Pi))). (A8) 

Here, q(t,qi,pi) and p(t,qi,pi) are the solution of the quasiclassical Langevin equation with the initial condition 
(q(0),p(0)) = (qi,Pi) and (} denotes the average over a(t). W n (qi,Pi,0) and W m (qi,Pi,T) are the initial and final 
Wigner functions. In order to utilize the result for the classical system <|A7|1 . we note that the Wigner function is very 

localized in semiclassical regime, (in fact for the harmonic potential case is proportional to e h mq where k 

is the strength of the potential). Therefore one may take a following view. Only the neighborhood of energy surfaces 
H(qi,pi, A(0)) = E n (0) and H(qf,pf,X(T)) = E m (T) does main contribution to the integration above. The work 
done on the system W for long enough time duration is nearly constant irrespective of the initial condition q(0),p(0). 
And the exponent of 1A7(I is within this approximation independent from the path a(t). Thus one can finally confirm 
the microscopic reversibility © for quasiclassical Langevin system. 

-Pf(KO) — >■ \m,T)) _ p( Em (T)-E n (p)-W+0(K)) ^ P -PQ (\q\ 
P R (e\m,T) ^ e\n,0)) 1 1 

Here the first equality results from the above view point and the O(H) deviation is caused by quasiclassical approx- 
imation due to the symmetrization procedure such as qp — > h(qp + pq) and this result is enough for the purpose of 
derivation of ©■ 
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